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NUMERICAL  MODELLING  OF  INTENSE  ELECTRON  BEAM  TRANSPORT 


IN  THE  SPIRAL  LINE  INDUCTION  ACCELERATOR 


I.  INTRODUCTION 

Advanced  high-current  accelerator  concepts,  such  as  the  heavy  ion 
accelerators  that  are  being  considered  as  inertial  fusion  drivers  [1],  or 
the  compact  high-current  electron  accelerators  [2]  that  are  under 
development  for  military  and  commercial  applications,  present  unique  and 
challenging  beam  transport  problems.  In  both  cases,  the  beam  currents  are 
in  the  multi-kiloamp  regime,  where  space  charge  and  self-magnetic  forces 
play  a  major  role.  In  the  heavy  ion  fusion  accelerator,  both  the 
longitudinal  and  transverse  emittance  are  strongly  constrained  by  the 
requirement  that  the  ion  beam  subsequently  be  focused,  over  a  focal  length 
of  several  meters  (the  radius  of  the  reaction  chamber),  onto  a  target  with 
a  radius  of  order  0.2  cm.  In  the  case  of  high-current  electron 
accelerators,  there  is  considerable  interest  in  generating  beams  of 
energetic  electrons  in  very  compact  devices.  To  minimize  size,  multiple- 
pass  designs  such  as  the  ion-focused  regime  recirculator[3] ,  the  modified 
betatron  accelerator  (MBA)[4]  and  the  spiral  line  induction  accelerator 
(SLIA) [5]  have  been  investigated.  Each  of  these  includes  tightly  curved 
transport  sections  with  radii  of  curvature  of  order  1  meter.  This  paper  is 
concerned  specifically  with  the  magnetic  transport  issues  that  arise  in  the 
context  of  SLIA. 

In  the  proposed  design  for  SLIA  (see  Fig.  1),  the  beam  is  transported 
along  an  open-ended  beam  pipe  which  makes  multiple  passes  through  each 
accelerating  cavity.  The  accelerating  cavities  are  aligned  along  the 
straight  sections  of  the  beam  pipe,  where  the  beam  is  transported  on  a 
solenoidal  focusing  field.  To  bring  the  beam  back  around  after  each 
straight  section,  it  is  taken  through  a  180°  achromatic  bend,  where  a 
strong  focusing  stellarator  (helical  quadrupole)  field  is  provided,  in 
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addition  to  the  solenoidal  field  and  the  vertical  bending  field.  Because 
the  field  coils  are  wound  directly  onto  the  spiral  beam  pipe,  and  because 
each  bend  is  magnetically  shielded  from  its  neighbors,  each  bend  can  be 
individually  designed  to  be  achromatic  near  the  expected  beam  energy  at 
that  point.  (A  perfectly  achromatic  bend  would  be  defined  as  one  in  which 
each  beam  particle  exits  the  bend  along  the  trajectory,  relative  to  the 
axis,  with  which  it  entered  the  bend.)  A  schematic  diagram  of  the  proposed 
SLIA  proof-of-concept  experiment  (P0CE)[5J,  which  uses  three  bends  and 
makes  two  passes  though  each  set  of  accelerating  cores,  is  shown  in  Fig.  1. 

Two  key  problems  must  be  addressed  in  this  transport  scheme.  The 
first  concerns  the  effect  of  the  hoop  forces  which  act  on  the  beam  in  the 
bends.  These  are  the  net  electrostatic  and  magnetostatic  forces  acting  on 
a  given  beam  segment,  due  to  the  remainder  of  the  curved  beam.  Hoop  forces 
have  been  analyzed  only  for  a  coasting  beam  in  a  circular  torus,  such  as 
the  MBA. [6]  The  second  difficulty  arises  from  the  fact  that  a  matched  beam 
in  a  solenoid  has  a  circular  cross-section,  while  a  matched  beam  in  the 
combined  solenoid/stellarator  fields  has  an  elliptical  cross-section. [7] 
Thus  it  is  necessary  to  use  matching  elements  at  the  entrance  and  at  the 
exit  of  each  bend,  in  order  to  avoid  mismatch  oscillations  of  the  beam 
cross-section  which  can  result  in  emittance  growth.  In  an  higher  energy 
SLIA  design,  there  can  be  on  the  order  of  10  passes  through  the 
accelerating  cores[2],  implying  19  bends  and  38  such  transitions. 

In  this  study,  we  focus  our  attention  on  two  specific  parameter 
regimes  of  immediate  interest:  (1)  the  POCE,  in  which  a  10  kA,  35  ns  beam 
is  to  be  injected  at  3.5  MeV  and  then  accelerated  tc  9.5  MeV  in  steps  of 
1.5  MeV,  and  (2)  recent  experiments[8-10]  in  which  an  850  keV,  200  A 
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injector  was  used  to  study  the  capture  of  the  beam  onto  a  solenoid,  and  the 
matching  at  the  entrance  of  the  stellarator  field. 

These  studies  are  based  on  numerical  simulation.  Initially,  the 
simulations  were  performed  using  the  ELBA  three-dimensional  particle 
simulation  code[ll],  which  solves  Maxwell's  equations  self-consistently 
with  the  complete  relativistic  beam  dynamics.  In  the  ELBA  code,  the  beam 
pipe  is  represented  as  a  straight  metallic  cylinder,  but  bends  are  included 
via  a  mapping  scheme  which  was  adapted  from  the  work  of  Friedman[12] ,  with 
corrections  to  the  electromagnetic  fields  which  are  valid  to  first  order  in 
the  small  parameter  ry/R.  Here,  rv  is  the  beam  pipe  radius  (3  cm  in  the 
POCE)  and  R  is  the  major  radius  of  the  bend  (80  cm  in  the  POCE).  The 
electromagnetic  field  corrections  are  discussed  in  the  Appendix.  The  bulk 
of  the  studies  reported  in  this  paper  were  done  with  a  faster  running 
version  of  the  code  called  SST  (Single  Slice  Transport) . 1 13]  Here,  self¬ 
field  effects  due  to  variations  in  the  beam  density  or  current  along  the 
length  of  the  beam  are  neglected,  so  that  only  a  single  transverse  slice  of 
the  beam  need  be  simulated.  For  transport  in  a  bend,  results  from  the  two 
codes  for  beam  displacement  and  optimum  vertical  field  showed  agreement  to 
within  less  than  IX  for  POCE  parameters. 

The  organization  of  the  paper  is  as  follows.  In  Section  II,  we 
discuss  the  initial  equilibrium  of  the  intense  beam  in  the  solenoid  field. 
In  section  III,  we  present  transport  simulations  for  each  of  the  three  POCE 
bends.  The  optimum  vertical  field  is  determined  for  each  bend  and  compared 
to  the  theoretical  estimate  of  Chernin  and  Sprangle.[6]  The  effect  of  this 
correction  to  the  vertical  field  on  the  achromatic  properties  of  the  first 
bend  is  also  assessed.  In  section  IV,  we  present  simulations  of  transport 
into  and  out  of  the  stellarator  fields  with  various  beam  profiles  for  POCE 
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parameters.  Ve  show  that  emittance  grovth  results  from  phase-mix  damping 
of  mismatch  oscillations  of  the  beam[ 14-16],  but  can  be  reduced  or 
eliminated  though  the  use  of  magnetic  focusing  elements.  Ve  present  a 
novel  method  of  accomplishing  this  match  via  a  single  thick  quadrupole 
lens.  Ve  present  end-to-end  simulations  of  the  POCE  with  and  without 
matching  elements.  In  section  V,  ve  discuss  the  recent  experimental 
results  of  Refs.  8-10.  Simulations  with  these  parameters  illustrate  the 
sensitivity  of  the  mismatch  oscillation  characteristics  to  the  radial 
distribution  of  the  beam.  Our  conclusions  are  given  in  section  VI. 
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II.  EQUILIBRIA  FOR  SOLENOIDAL  TRANSPORT 

In  accelerator  transport  studies,  it  is  frequently  assumed  that  the 

beam  has  a  Kapchinskij -Vladimirskij  (K-V)  distribution! 17]  in  the  rotating 

frame  in  the  directions  (x,y)  transverse  to  the  propagation  direction  (z). 

The  K-V  distribution,  which  has  a  flat  radial  density  profile,  is  an  exact 

equilibrium  in  the  limit  y^  •+  ®  and  v/y^  ■*  0,  where  y^  is  the  relativistic 

factor  corresponding  to  the  beam  kinetic  energy,  v  =  I^[kA]/17B^  is 

2  1/2 

Budker's  parameter,  1^  is  the  beam  current  and  0^  ■  (1  -  1/y^) 

However,  particle  simulations  of  the  POCE  beam  in  the  5.492  kG  solenoidal 

guide  field  show  ar  instability  if  the  beam  is  injected  with  a  transverse 

K-V  distribution.  The  beam  radius  exhibits  growing  oscillations  at  the 

envelope  frequency  which  typically  reach  amplitudes  of  the  order  of  7  to 

102  of  the  equilibrium  radius.  These  are  shown  in  Fig.  2  (solid  line), 

where  the  RMS  radius  r^  is  plotted  versus  transport  distance  z  for  the  POCE 

beam  (10  kA,  4  MeV,  normalized  RMS  emittance  e  «  0.158  cm-rad) 

n,  rms 

initialized  with  a  K-V  distribution.  Note  that  in  this  simulation  and 
throughout  the  paper  we  assume  that  the  diode  is  nonimmersed,  so  the 
canonical  angular  momentum  of  the  beam  is  zero.  Eventually  (over  a  scale 
length  of  several  meters),  the  oscillation  amplitude  saturates.  The 

instability  is  a  kinetic  effect  which  does  not  occur  in  envelope  code 
calculations  and  is  not  suppressed  by  simply  initiating  the  beam  with 
carefully  matched  envelope  parameters.  It  indicates  that  the  initial 

distribution  is  not  a  true  equilibrium,  and  that  the  beam  dynamics  are 

acting  to  rearrange  the  distribution  to  approach  equilibrium.  Similar 
phenomena  are  seen  if  the  beam  is  initialized  in  other  distributions  that 
are  not  exact  equilibria  for  finite  v/y^,  e.g.,  a  parabolic  density  profile 
with  a  Gaussian  velocity  distribution.  Although  the  instability  is  rather 
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weak,  it  tends  to  mask  the  features  of  interest.  Oscillations  of  this  type 
are  probably  inevitable  in  a  real  beam,  where  one  does  not  have  perfect 
control  over  the  beam  distribution. 

The  K-V  distribution  can  be  derived  from  the  Vlasov  equilibrium  given 

2 

by  f(r,p)  ■  nQ/2nrQm  $(H-wL-YqIiic  )  SCP-BqYqIhc)  in  the  limit  Yq  -*  ®.  [18] 
Here  the  energy  H,  the  canonical  angular  momentum  L,  and  the  canonical 
axial  momentum  P  are  constants  of  motion  and  »,  Yq»  3q  and  n^  are 
parameters.  We  have  derived  an  improved  distribution  by  retaining  the 
first  order  terms  in  1/Yq»  which  result  in  corrections  to  the  K-V  including 
the  self-generated  diamagnetic  axial  magnetic  field  and  the  kinetic  energy 
variations  due  to  the  space  charge  depression. [ 19]  Although  the  charge 
distribution  is  not  radially  flat  in  this  approximation,  it  has  a  sharp 
cutoff  as  does  the  K-V.  When  the  POCE  beam  is  initialized  with  this 
higher-order  K-V  distribution,  as  in  Fig.  2  (dashed  line),  growing 
oscillations  of  the  beam  radius  are  not  observed. 

Simulations  of  the  POCE  injector[20] ,  as  well  as  experimental 
measurements[8-10] ,  indicate  that  the  beam  injected  into  the  SLIA  transport 
system  has  a  rounded  radial  profile,  and  thus  does  not  have  a  K-V  type 
distribution.  In  simulation  studies,  we  have  found  it  convenient  to  use  a 
numerical  technique  to  initiate  the  beam  equilibrium,  rather  than  the 
analytic  expansion  described  in  the  previous  paragraph.  We  begin  the 
simulation  with  a  beam  distribution  that  reflects  the  macroscopic 
properties  of  interest  (radial  density  profile  and  emittance),  and  with 
carefully  matched  envelope  parameters,  including  an  ad  hoc  correction  to 
the  initial  particle  energy  as  a  function  of  radius  due  to  the  space  charge 
potential  of  the  beam: 

\>.  2  2.  .  2  ... 

Y  ”  Yb  +  2(r  "  rb)/rb  ’  (1) 
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where  y  is  the  relativistic  factor  assigned  to  a  particle  at  radius  r  and 
r^  is  the  RMS  radius  of  the  beam.  We  have  followed  this  procedure  for  a 
variety  of  equilibria  and  find  that  the  space-charge  correction  reduces  but 
does  not  eliminate  the  unstable  oscillations.  The  beam  is  then  transported 
numerically  for  several  meters  in  the  solenoid  field,  thus  allowing  the 
oscillations  to  saturate  and  damp.  Such  equilibria  were  used  to  initiate 
the  simulations  presented  below. 
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III.  TRANSPORT  THROUGH  THE  180°  BEND 


Hoop  forces  occur  in  the  bends  because  a  given  point  in  the  curved 
beam  experiences  a  net  transverse  electro-  and  magnetostatic  force  due  to 
the  remainder  of  the  curved  beam.  These  self-fields  have  been  analyzed  in 
some  detail  for  a  beam  of  circular  cross-section  in  a  toroidal 
configuration. [6]  Assuming  the  beam  is  on  axis,  the  optimum  vertical  field 
is  given  by 


2 

mrb*bc  1  me2  \>l\,  ,„21 

eR -  +  \~r  rL(1  +  eb)ln(?T)  +  **bj 


(2) 


where  m  is  the  electron  mass,  R  is  the  major  radius  of  the  bend,  ry  is  the 
radius  of  the  beam  pipe,  e  (taken  to  be  positive)  is  the  elementary  charge 
and  c  is  the  speed  of  light.  The  parameter  6  depends  on  the  details  of  the 
beam  equilibrium  in  the  bend:  6  =  1  if  «ne  assumes  consta  .t  current 

density  across  the  beam  cross-section  and  6  *  2  if  one  assumes  constant 

circulation  frequency.  The  first  term  in  Eq.  (2)  is  the  ideal  value  for  a 
low-current  beam;  the  second  is  required  to  correct  for  the  hoop  forces. 

For  the  POCE,  in  which  I,  *  10  kA,  R  =  80  cm  and  r  =  3  cm,  these 

D  V 

forces  can  be  substantial.  In  the  first  bend,  for  example,  the  design 
energy  is  4  MeV  and  the  ideal  vertical  field  is  186.9  G.  Equation  (2) 
gives  Bv  between  236.2  G  (for  6*1)  and  248.6  G  (for  6=2).  There  are 
additional  uncertainties  in  these  values  because  a)  the  beam  cross-section 
in  the  bend  is  not  circular  but  is  instead  a  slowly  rotating  ellipse  and  b) 
the  SLIA  bend  extends  for  180®  with  entrance  and  exit  points,  whereas  the 
published  calculation  applies  to  a  full  toroidal  configuration. 

Field  parameters  for  the  three  POCE  achromatic  bends  are  given  in 
Table  1.(21)  In  Table  I,  Bg  i:  the  axial  magnetic  field  at  the  geometric 
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axis,  bq  is  the  quadrupole  field  gradient  near  the  quadrupole  field  axis, 
and  Xp  is  the  pitch  length  of  the  quadrupole  windings.  The  SST  and  ELBA 
codes  use  the  following  somewhat  idealized  approximations  for  the  external 
fields  in  the  bend.  The  axial  field  is  given  by 


B  «  B  R/(R  +  x) 
z  s 


The  quadrupole  fields  are 


B  -  b  [y  cos(k  z)  -  (x  -  x  )sin(k  z)] 
a  q  q  q  q 


By  *  bq[(x  -  xq)cos(k^z)  +  y  sin(kqz)]  , 


where  k  *  4n/X  and  x  is  the  x-coordinate  of  the  quadrupole  field  axis, 
q  p  q 

which  is  displaced  relative  to  the  geometric  axis  in  the  bend[22]: 


■ft  *  £-f] 


The  vertical  field  enters  as 


By  -  -Bv(l  -  rpc/R) 


Bx  *  Bvr*y/,R  * 


where  h  is  the  vertical  field  index  and  is  responsible  for  weak  focusing  in 
betatrons.  In  this  coordinate  system,  the  beam  moves  in  the  z  direction, 
and  x  and  y  are  the  horizontal  and  vertical  coordinates,  respectively.  The 
curve  lies  in  the  x-z  plane  with  the  inside  wall  of  the  curved  beam  pipe  at 


x  =  -r  .  In  Eqs.  (3-6),  we  assume  rv/R  \<  1. 
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Each  bend  is  designed  to  be  achromatic  near  the  expected  beam  energy 

for  that  bend.  In  other  words,  if  the  beam  is  within  10-20%  of  the  design 

energy  for  a  given  bend,  the  beam  should  exit  with  the  same  position  and 

angle,  relative  to  the  axis,  that  it  had  when  it  entered.  A  perfectly 

achromatic  bend  will  have  no  effect  on  the  beam  trajectory  relative  to  the 

axis  no  matter  what  the  beam  energy.  In  practice,  parameters  have  been 

determined  that  are  correct  to  first  order  in  the  small  parameter  AE/E 

where  E  is  the  design  energy  for  the  bend  and  AE  is  the  deviation  of  the 

beam  energy  from  the  design  energy. [21,23] 

Simulation  results  for  the  first  bend  show  that  if  Bv  is  set  to  the 

value  By  «  186.9  G  which  includes  no  hoop  force  corrections,  the  beam 

centroid  displacement  within  the  bend  will  exceed  1  cm.  This  is  shown  in 

Fig.  3  where  the  x  and  y  positions  of  the  beam  centroid  are  plotted  versus 

transport  distance  for  the  first  POCE  bend.  Here  the  bend  entrance  is 

located  at  z  «  50  cm  with  the  exit  at  z  *  301.3  cm.  Even  with  this 

unacceptably  large  displacement,  the  achromatic  design  brings  the  centroid 

back,  to  within  2  mm  of  the  axis  on  exit.  A  series  of  runs  were  performed 

to  determine  the  maximum  beam  centroid  displacement  r  as  a  function  of 

c  j  max 

the  imposed  Bv-  The  results  are  shown  in  Fig.  4.  The  optimal  value  of 
the  vertical  field  is  seen  to  be  233  G,  for  which  the  beam  centroid  remains 
within  0.06  cm  of  the  geometric  axis  at  all  times.  This  value  of  B^  is  a 
25%  increase  over  the  ideal  value  but  is  somewhat  smaller  than  expected 
from  Eq.  (2).  Numerically  determined  optimum  vertical  field  values  for  all 
three  bends  are  given  in  Table  II.  For  comparison,  the  ideal  values  (no 
hoop  force  correction)  and  the  predicted  values  from  Eq.  (2)  are  also 
given. 
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With  a  25%  increase  in  the  design  value  of  the  vertical  field,  it 
might  be  expected  that  other  corrections  to  the  design  parameters  of  the 
stellarator  and  solenoidal  fields  would  be  necessary  to  preserve  the 
achromatic  property  of  the  bend.  For  instance,  an  achromatic  design 
requires  that  the  frequencies  associated  with  the  various  fields  produce 
integer  numbers  of  oscillations  as  the  beam  traverses  the  bend.  However, 
Mondelli  et  al.[24]  have  shown  that  the  vertical  field  correction,  with  an 
appropriate  correction  to  the  vertical  field  index,  simply  preserves  the 
frequency  that  would  be  expected  in  the  ideal,  low-current  case.  Figures  5 
and  6  confirm  this  conclusion.  These  figures  show  the  centroid  position 
xc(z),  yc(z)  as  the  beam  passes  around  the  first  POCE  bend  with  the  optimum 
vertical  field  value  given  in  Table  I.  In  Fig.  5,  the  beam  energy  was 
offset  from  the  design  energy  by  AE  *  -10%.  In  Fig.  6,  AE  «  +10%.  In  both 
cases,  the  beam  centroid  returns  to  within  0.2  cm  of  the  axis. 

In  principle,  it  is  also  necessary  to  correct  the  vertical  field  index 
q  to  account  for  diamagnetic  effects  in  the  intense  beam,  which  tend  to 
reduce  h  to  less  than  the  design  value  over  most  of  the  beam  cross-section. 
We  have  performed  several  simulations  with  h  set  between  0  and  1.0  (at  low 
current,  the  design  value  is  t)  ■  0.5).  These  simulations  showed  that  the 
value  of  the  vertical  field  index  had  less  than  a  1%  effect  on  the  beam 
centroid  position,  and  had  no  effect  whatsoever  on  transport  at  the  design 
energy  or  the  achromatic  properties  of  the  first  bend.  This  result  might 
well  have  been  anticipated  inasmuch  as  the  strong  focusing  due  to  the 
stellarator  fields  dominates  the  weak  focusing  arising  from  the  betatron 
field.  In  subsequent  simulations,  we  simply  set  h  =  0.5. 

We  also  measured  the  sensitivity  of  the  achromatic  design  to  errors  in 
the  settings  of  the  axial  and  stellarator  magnetic  fields.  In  Figs.  7  and 
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8,  we  show  the  results  of  setting  AE  *  -10%,  and  varying  B  and  b  .  It  is 

S  q 

apparent  that  the  achromatic  design  is  tolerant  to  errors  in  the  2-3%  range 
and  that  the  design  values  are  not  optimized  for  AE  =  -10%.  Further  runs 
at  AE  <  + 10 X  showed  similar  sensitivity  to  field  errors.  The  optimum 
values  at  AE  ®  +10%  differ  from  those  at  AE  =  0  or  AE  *  -10%  but  not  in  an 
easily  predictable  way  (i.e.,  corrections  to  the  optimum  fields  at  AE  = 


+10%  are  not  simply 

the  negatives 

of  those 

at  AE  = 

-10%). 

These 

observations  are  not 

surprising  in 

that  the 

achromatic 

parameters  are 

correct  only  to  first 

order  in  AE/E. 
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IV.  MATCHING  VIA  MAGNETIC  FOCUSING  ELEMENTS 

A  major  concern  in  the  SLIA  design  is  the  transition  at  the  entrance 
and  exit  of  each  bend.  In  the  straight  sections,  the  beam  is  transported 
on  a  solenoidal  field  in  an  equilibrium  state  with  a  circular  cross- 
section.  In  the  bends,  a  strong  focusing  stellarator  (rotating  quadrupole) 
field  is  also  imposed,  to  improve  momentum  compaction.  Chernin[25]  has 
derived  an  envelope  equation  for  transport  in  stellarator/solenoidal  fields 
and  has  shown  that  matched  solutions  exist  in  which  the  beam  cross-section 
is  elliptic,  with  constant  major  and  minor  radii  and  with  the  orientation 
of  the  beam  ellipse  rotating  in  phase  with  the  stellarator  windings.  Since 
the  transitions  onto  and  off  of  the  bends  will  be  made  more  than  fifty 
times  in  a  50  to  100  MeV  accelerator,  it  is  crucial  that  each  transition  be 
made  smoothly  and  with  a  minimum  of  emittance  growth. 

For  beams  with  coupled  x-y  motion,  Chernin[25]  has  provided  a  suitable 
definition  of  emittance  which  is  a  straight-forward  generalization  of  the 
usual  emittance  definition  (e.g. ,  Ref.  26).  This  definition  is 


e 


rms 


(7) 


where  £  is  a  4  by  4  matrix  with  elements 

E.j  «=  <u.Uj>  -  <u.><Uj>  ,  (8) 

Uj  ■  x,  U2  *  dx/dz,  u^  ■  y>  u4  *  dy/dz  and  "<  >"  signifies  an  average  over 

the  beam  particles.  The  emittance  t  is  conserved  by  Chernin's  envelope 

rms 

equation[25] ,  which  makes  the  usual  assumption  that  all  forces  acting  on 
the  beam  are  linear  in  the  radial  coordinate  r.  It  is  therefore  assumed 
that  the  beam  density  profile  is  flat-topped,  which  is  necessary  if  the 
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space-charge  forces  are  to  be  linear.  The  normalized,  RMS  emittance  is 

e  _  -  $yt  .  In  terms  of  the  matrix  elements  I..,  the  major  radius  r 

n,  rms  rms  ij  + 

and  the  minor  radius  r  of  the  beam  cross-section  are  given  by 


Z11  +  Z33  *  r0 


(9) 


2  2  1/2 

where  Iq  «  ((E^  -  £33)  +  ^E^]  '  For  a  beam  yith  circular  cross- 

section,  r  *  r  *  4T  r,  . 

+  -  D 

Simulations  show  that  large-amplitude  envelope  oscillations  are 
excited  at  the  entrance  and  exit  of  the  first  POCE  bend  if  no  matching  lens 
elements  are  used  to  facilitate  the  transition.  The  characteristics  of  the 
oscillations  (wavelength,  damping  and  emittance  growth)  are  not  determined 
entirely  by  beam  envelope  parameters  but  are  sensitive  to  the  distribution 
function  of  the  beam. 

In  Fig.  9,  we  show  r^(z)  for  the  POCE  beam  as  it  transits  the  first 
bend  of  the  POCE.  Four  different  beam  profiles  are  simulated.  In  Fig. 
9(a),  the  beam  is  initiated  with  a  K-V  distribution.  At  the  entrance  to 
the  bend  (z  «  50  cm),  the  beam  transitions  to  an  elliptic  cross-section, 
but  overshoots  the  elliptic  equilibrium  and  enters  into  a  series  of 
undamped  large-amplitude  oscillations,  which  are  particularly  evident  in 
r+(z).  At  the  exit  from  the  bend  (z  =  301.3  cm),  the  beam  transitions  to 
an  elliptic  cross-section  of  lower  eccentricity.  After  a  series  of  damped 
oscillations  of  r+(z),  the  beam  approaches  an  elliptic  equilibrium  in  the 
solenoidal  transport  system,  but  shows  no  sign  of  damping  back  toward  the 
conventional  equilibrium  with  circular  cross-section.  Little  or  no 
emittance  increase  occurs,  as  seen  in  Fig.  10. 

This  type  of  behavior  is  characteristic  of  the  K-V  distribution. 
Because  the  radial  profiles  of  charge  and  current  density  are  nearly  flat, 
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the  self-forces  are  nearly  linear,  there  is  very  little  phase-mix  damping, 
and  Chernin's  envelope  equations[25]  provide  a  nearly  complete  description 
of  the  beam  dynamics.  The  damping  of  the  weak  oscillations  about  an 
elliptic  equilibrium  in  the  solenoid  is  probably  associated  with  self-force 
nonlinearities  that  result  from  the  inexactness  of  the  K-V  equilibrium,  as 
discussed  in  Sec.  II. 

The  results  are  markedly  different  if  the  beam  profile  is  rounded,  as 
can  be  seen  in  Figs.  9(b)-9(d).  In  case  9(b),  the  beam  is  initiated  with  a 
spatially  parabolic  profile, 


nb<r> 


'  nb(r.0)[l  -  <r/2rb)2]  , 
.  0 


r  *  2rb 


r  >  2r 


b  ’ 


(10) 


where  rfa  is  the  RMS  beam  radius.  In  Fig.  9(c),  the  beam  is  initiated  with 
a  Gaussian  spatial  profile  and,  in  Fig  9(d),  with  a  truncated  Bennett[27] 
spatial  profile, 


nb<r> 


f  nb(r=0)[l  +  (r/r2)}  2 


l  0 


r  <  Srh/JT"1 
r  >  3rb/'f2_l  . 


(ID 


Each  of  our  non-K-V  beams  have  Gaussian  distribution  functions  in  velocity 
space.  In  each  of  these  cases,  the  beam  becomes  mismatched  at  the 
stellarator  entrance,  damps  toward  an  elliptical  state  in  the  bend,  becomes 
mismatched  again  at  the  exit,  and  finally  damps  toward  a  circular  cross- 
section.  Each  stage  of  damping  is  accompanied  by  emittance  growth  as  seen 
in  Fig.  10.  The  damping  would  appear  to  be  associated  with  phase  mixing, 
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due  to  the  nonlinear  self-field  forces  arising  from  the  rounded  beam  radial 
profiles.  Since  rounded  beam  profiles  are  typical  in  the  SLIA,  the 
pictures  presented  in  Figs.  9(b-d)  are  representative  of  SLIA  physics, 
rather  than  the  K-V  case,  Fig  9(a),  where  the  oscillations  are  highly 
coherent  and  undamped.  Ve  can  expect  considerable  emittance  growth  in  the 
accelerator  unless  lens  elements  are  introduced  to  carefully  match  the  beam 
at  the  entrance  and  exit  of  each  bend. 

Tiefenback[9]  has  shown  that  the  beam  envelope  parameters  can  be 
matched  exactly  by  using  a  compound  magnetic  lens  consisting  of  three  ideal 
thin  focusing  elements  (one  solenoid  lens  and  two  quadrupole  lenses). 
Petillo  et  al . [ 21 )  have  used  an  envelope  code  to  show  that  this  scheme  can 
be  implemented  with  realistic,  finite  thickness  lenses  to  obtain  a  useful 
inexact  match.  Ve  have  found  a  novel  solution  to  the  matching  problem 
which  employs  a  single  thick  quadrupole  lens. 

The  key  to  this  matching  scheme  is  the  effect  of  a  thick  quadrupole 
lens  on  beam  transport  in  a  solenoidal  field.  This  is  illustrated  in  Fig. 
11,  where  we  plot  the  major  and  minor  radii  r±(z)  obtained  from  solutions 
of  the  envelope  equations  with  beam  and  solenoidal  field  parameters 
corresponding  to  those  of  the  first  POCE  bend.  A  thin  quadrupole  lens, 
Fig.  11(a),  simply  induces  a  series  of  quasiperiodic  quadrupole  envelope 
oscillations,  i.e.,  the  beam  cross-section  oscillates  between  states  of 
maximum  elliptic  eccentricity,  passing  through  a  circular  cross-section  in 
between  maxima.  However,  when  the  quadrupole  lens  has  finite  thickness  Lq 
[defined  as  the  full-width-at-half-maximum  length  of  the  field  envelope 
bQ(z)],  as  in  Fig.  11(b),  the  traces  of  r+(z)  and  r_(z)  separate,  so  that 
the  cross-section  is  non-circular  at  all  times,  but  oscillates  between 
states  of  maximum  and  minimum  elliptic  eccentricity.  As  the  thickness  of 
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the  lens  increases  further,  to  Lq  -  5  cm  in  Fig.  11(c),  r+(z)  and  r  (z) 
become  veil  separated,  and  oscillate  separately  around  central  values 
corresponding  to  an  elliptic  equilibrium  state.  For  a  still  thicker  lens, 
Lq  •  10  cm  in  Fig.  11(d),  r+(z)  performs  small-amplitude  sinusoidal 
oscillations  about  a  well-defined  equilibrium  value  and  phase  locking  is 
observed  such  that  r  (z)  is  approximately  stationary  at  the  extrema  of 
r+(z).  The  beam  has  been  transformed  from  its  initial  circular  equilibrium 
to  a  new  approximate  equilibrium  in  which  the  cross-section  is  elliptic  and 
slowly  rotating.  The  eccentricity  of  the  ellipse  is  nearly  constant  and  is 
determined  by  the  integrated  strength  of  the  lens  Jdz  bQ(z).  Recall, 
however,  that  in  the  context  of  a  more  exact  kinetic  treatment,  rather  than 
an  envelope  model,  the  elliptic  equilibrium  will  eventually  damp  back  to 
circular  if  there  is  phase  mixing  due  to  a  rounded  beam  radial  profile  (see 
Fig.  9).  The  transition  to  an  elliptic  equilibrium  via  a  single  quadrupole 
lens  and  subsequent  damping  toward  a  circular  beam  cross-section  has  been 
observed  both  experimentally[22]  and  in  particle  simulations. 

Envelope  modelling  and  particle  simulations  indicate  that  the  beam 
response  to  a  quadrupole  lens  is  not  very  sensitive  to  the  details  of  the 
quadrupole  lens  axial  dependence  bQ(z),  but  rather  depends  primarily  on  the 
integrated  lens  strength  Jdz  bQ(z)  and  on  the  thickness  Lq.  An  approximate 
criterion  for  the  "thick  lens"  limit,  as  exemplified  by  Fig.  11(c)  and  to 
an  even  greater  degree  by  Fig.  11(d),  is 


X  „  2 

L  >  _c  . 

LQ  4  2eB 


(12) 


where  X  is  the  cyclotron  wavelength  in  the  solenoidal  field.  (For  SLIA 
c 

parameters,  the  oscillation  wavelength  is  not  significantly  increased  by 
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space-charge  fields.)  In  the  first  POCE  bend,  X£  *  17.2  cm,  and  the  lens 
is  "thick"  in  Figs.  11(c)  and  11(d),  vhere  Lq  -  5  cm  and  10  cm, 
respectively. 

The  matching  onto  the  stellarator  equilibrium  is  accomplished  by  using 
a  quadrupole  lens  that  is  in  the  "thick  lens"  limit  (given  the  beam  energy 
at  that  point  in  the  beam  line),  adjusting  the  integrated  strength  of  the 
lens  to  match  the  eccentricity  needed  in  the  bend,  adjusting  the  distance 
between  the  lens  and  the  bend  such  that  the  beam  enters  the  stellarator 
field  at  a  stationary  point  of  the  oscillation,  and  adjusting  the 
orientation  of  the  lens  so  that  the  angle  of  the  ellipse  is  in  the  correct 
phase  as  the  beam  enters  the  bend.  A  single  quadrupole  lens  can  also  be 
used  at  the  exit  of  the  bend  to  bring  the  beam  back  to  a  nearly  circular 
cross-section.  In  the  ideal  case  (no  envelope  oscillations),  the  exit 
quadrupole  is  identical  to  the  entrance  quadrupole. 

It  should  be  noted  that  the  required  minimum  lens  thickness  Lq,  such 
that  the  lens  is  "thick,"  increases  proportionally  to  the  beam  energy, 
according  to  Eq.  (12).  For  example,  the  lens  of  Fig.  11(c)  is  sufficiently 
thick  for  the  first  bend  of  the  POCE  (electron  energy  4  MeV),  but  not  for 
the  third  bend  (7  MeV),  as  shown  in  Fig.  12.  For  practical  reasons,  it  may 
be  necessary  to  use  another  matching  scheme  at  high  energy. 

Although  the  features  of  the  matching  scheme  are  easily  accessed  and 
surveyed  by  numerically  solving  the  envelope  equations,  we  have  used 
particle  simulation  to  study  the  transport  processes  in  full  detail. 
Figures  13  and  14  show  the  results  of  a  simulation  of  transport  around  the 
POCE  first  bend  for  a  case  in  which  the  beam  is  initiated  with  a  parabolic 
radial  profile  identical  to  that  of  Fig.  9  (b).  Here,  a  single  quadrupole 
lens  is  used  to  match  the  beam  at  the  entrance  to  the  bend.  The  lens  is 
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located  25  cm  from  the  stellarator  entrance,  with  peak,  quadrupole  gradient 
bg  -  292  G/cm,  z-dependence  as  plotted  in  Fig.  15  and  lens  thickness  Lq  - 
5.5  cm.  At  the  bend  exit,  ve  fine-tuned  the  match  using  an  exit  quadrupole 
with  bg  *  255  G/cm  (this  fine  tuning  accounts  for  the  phase  of  the  small 
mismatch  oscillation  as  the  beam  exits  the  stellarator).  The  lens  profile 
of  Fig.  15  is  taken  from  measurements  of  an  actual  quadrupole  lens  used  in 
SLIA  experiments. [22]  Vith  this  matching  scheme,  the  radial  mismatch 
oscillations  of  Fig.  9(b)  have  been  minimized  and  the  13%  emittance  growth 
of  the  unmatched  case  has  been  eliminated.  For  comparison,  Fig.  14  shows 
emittance  plotted  versus  transport  distance  for  this  case  and  for  the 
simulation  of  Fig.  9(b).  We  have  obtained  similar  results  via  the  three- 
element  match  of  Ref.  21. 

These  simulations  have  been  extended  to  provide  complete  end-to-end 

particle  simulations  of  transport  in  the  POCE  with  each  of  the  three  bends 

modeled  as  described  above.  The  acceleration  sections  are  modeled  by 

simply  imposing  an  axial  electric  field  over  a  short  section  of  the  beam 

pipe.  The  location  of  each  curve  and  each  idealized  acceleration  section 

for  the  end-to-end  simulations  is  given  in  Table  III.  Results  from  a 

simulation  without  matching  elements  are  shown  in  Fig.  16.  The  beam  is 

initiated  with  a  parabolic  spatial  profile  as  in  Figs.  9(b)  and  13  above. 

Overall,  the  emittance  t  increases  by  52%.  Plots  of  r  (z)  show  that 

n  f rms  t 

the  beam  becomes  mismatched  in  each  bend  and  damps  back  to  an  approximately 
circular  cross-section  after  each  bend.  It  is  interesting  to  note  that,  by 
chance,  the  beam  exits  the  third  bend  with  an  approximately  circular  cross- 
section  and  little  subsequent  emittance  growth  occurs.  For  comparison,  we 
repeated  this  simulation  with  the  third  bend  terminated  at  z  =  1250  cm,  at 
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the  point  of  maximum  eccentricity  of  the  beam  cross-section.  In  this  case 
ve  observed  an  overall  emittance  increase  of  110%. 

The  results  of  an  end-to-end  POCE  simulation  in  which  single-element 
matching  is  employed  on  all  three  bends  are  shown  in  Fig.  17.  Here,  the 
initial  beam  parameters  are  identical  to  those  of  Fig.  16.  With  the  use  of 
single-element  matching,  however,  the  normalized  emittance  is  largely 
preserved,  increasing  by  only  10%.  In  this  simulation,  quadrupole  lenses 
with  peak  quadrupole  gradient  b^  *  149  G/cm  were  placed  4  cm  from  the 
entrance  and  4  cm  after  the  exit  of  the  second  bend.  For  the  third  bend, 
lenses  were  located  6  cm  from  the  stellarator  end  points  with  b^  «  266 
G/cm.  The  quadrupole  lens  profile  of  Fig.  15  was  used  for  the  first  bend, 
as  in  Fig.  13,  and  a  similarly  shaped  quadrupole  with  Lq  ■  11  cm  was  used 
on  the  second  and  third  bends. 
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V. 


SIMULATION  OF  THE  850  KeV  EXPERIMENT 


In  preparation  for  the  POCE,  a  series  of  preliminary  transport  and 

matching  experiments  has  recently  been  performed  at  PSI.[8-10]  The 

parameters  for  these  experiments  were  chosen  to  replicate  as  closely  as 

possible  the  physics  of  the  POCE  beam  at  much  lower  current  and  energy. 

2  2 

Thus,  the  dimensionless  ratio  a  =  vr.  /yt  ,  which  determines  the  relative 

b  n 

influences  of  the  emittance  and  space  charge  forces,  was  chosen  to  be  the 
same  for  both  beams.  The  significance  of  a  can  be  noted  from  the  envelope 
equation  for  a  beam  in  the  presence  of  an  axial  magnetic  field[28,15] , 


2  2  2  b 
c  rK 


2  v 


2  2  3 


b  b 


0  , 


(13) 


where  Q z  *  elWmc  is  the  nonrelativistic  cyclotron  frequency  and  en  = 

4e  .  The  nominal  POCE  parameters  are  I,  •  10  kA,  y.  =  8.83  (4  MeV), 

n , rms  b  d 

e  0.159  cm-rad  and  r.  «  0.5  cm,  which  gives  a  *  0.041.  The  low 

n, rms  b  ° 

energy  experiments  were  done  at  1^  *  200  A,  y^  *  2.66  (850  keV),  and  r^  = 
0.25  cm.  Preservation  of  the  ratio  a  requires  e  *  0.022  cm-rad  (the 

actual  experimental  value  was  somewhat  lower,  as  will  be  discussed  below). 
In  these  experiments,  the  beam  was  transported  along  a  1.5  kG  solenoid 
field  into  a  stellarator  section  with  length  251.3  cm  and  pitch  *  68.33 
cm,  identical  to  that  of  the  first  POCE  bend.  This  stellarator  section, 
however,  was  straight  (rather  than  curved)  so  as  to  facilitate  diagnostics, 
and  the  stellarator  strength  was  reduced  to  b  =68.5  from  the  POCE  value 
of  251.6  G/cm. 

In  the  initial  experiment,  the  beam  was  transported  into  the 
stellarator  fields  without  matching  at  the  stellarator  entrance.  The  data 
from  this  experiment,  giving  the  beam  major  and  minor  radii  r^(z),  are 
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shown  in  Fig.  18.  Each  data  point  represents  a  time-integrated  average 
over  one  shot  of  the  device,  obtained  by  inserting  a  piece  of  film  into  the 
beam  pipe  at  a  location  that  was  varied  from  shot  to  shot.  It  should  be 
noted  that  the  experimental  "major  and  minor  radii"  r+(z)  are  determined  by 
fitting  a  constant  current-density  contour  at  approximately  10%  of  the  peak 
current  density  in  each  shot,  whereas  the  simulation  r+(z)  values  are 
computed  via  Eq.  (9)  and  are  RMS  values  multiplied  by  'f2n.  The 
experimental  measurements  extend  only  to  z  «  350  cm,  50  cm  beyond  the  end 
of  the  stellarator  section.  As  a  result,  the  damping  of  mismatch 
oscillations  at  the  stellarator  exit,  as  well  as  the  emittance  growth 
associated  with  this  damping,  was  not  measured.  The  simulations  extend  to 
z  ■  500  cm,  enabling  us  to  study  matching  phenomena  at  the  ex_ t  as  well  as 
the  entrance  of  the  stellarator  section. 

In  Fig.  18,  it  is  evident  that  envelope  oscillations  are  excited  at 
the  stellarator  entrance  which  subsequently  damp  toward  an  elliptic 
equilibrium.  Emittance  measurements,  which  tend  to  be  problematical  for  a 
space-charge  dominated  beam,  were  made  both  at  the  injector  and  at  the  exit 
of  the  stellarator  section.  The  inherent  difficulties  associated  with  this 
measurement  were  further  exacerbated  at  the  and  of  the  stellarator  section 
by  the  need  to  aperture  the  beam  and  to  transport  the  resulting  "beamlet" 
into  a  field-free  region.  Because  the  measurement  was  performed  within  the 
3  cm  radius  beam  pipe,  the  drift  region  over  which  the  emittance  could  be 
measured  was  limited  to  10  cm.  In  fact,  a  definitive  measurement  of  the 
emittance  at  the  stellarator  exit  has  not  yet  been  obtained . [22 J 

Preliminary  measurements  of  the  beam  initial  radial  profile  showed  an 
approximately  parabolic  profile  with  emittance  en  rmg  =  0.008  rad-cm  and 
beam  radius  0.25  cm. [22]  In  our  first  simulation  studies,  we  initiated  the 
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beam  with  a  parabolic  profile.  One  such  simulation  is  shown  in  Fig.  19, 

with  beam  parameters  r,  «  0.25  cm,  e  *  0.011  rad-cm.  The  damping  rate 

b  n, rms 

was  found  to  be  comparable  to  that  of  the  experiment,  but  the  mismatch 
oscillation  frequency  differs  markedly  (by  252)  from  that  seen  in  the 
experiment.  It  was  found  to  be  impossible  to  match  the  experimental  data 
in  simulations  initiated  with  any  of  the  usual  distributions  (KV, 
parabolic,  Gaussian,  Bennett). 

Since  it  is  known[29]  that  core-halo  radial  profiles  are  associated 
with  increased  phase-mix  damping  of  radial  mismatch  oscillations  expected 
to  be  present  in  the  injector! 10, 20] ,  we  then  tried  initiating  the  beam 
with  a  core-halo  profile  in  the  simulations.  It  was  found  that  the 
simulations  could  provide  a  good  match  to  the  data.  Subsequently,  careful 
measurements  were  performed  which  indicated  that  the  beam  had  approximately 
102  of  its  current  in  a  diffuse  halo  extending  to  about  twice  the  radius  of 
the  main  core,  as  shown  in  Fig.  20. 

Creating  a  core-halo  distribution  for  simulation  purposes  was 
accomplished  by  adding  a  constant  density  spatial  distribution  with  edge 
radius  2r^  to  the  core  parabolic  distribution.  The  result  was  then 
transported  for  30  meters  and  allowed  to  come  to  equilibrium.  The  spatial 
profile  of  our  numerical  core-halo  distribution  is  shown  in  Fig.  21.  Here, 
we  put  252  of  the  current  into  the  initial  flat-top  halo.  In  the  resulting 
beam,  particles  with  r  >  r+  =  3.46  mm  represented  82  of  the  total  current. 
In  Fig.  21,  the  beam  emittance  is  t  *=  0.01  cm-rad.  If  we  consider 

only  the  particles  within  the  r  <  r  we  have  E  =  0.0065  cm-rad. 

Figure  22  shows  the  major  and  minor  radii  r+(z)  and  the  emittance 

t  (z)  of  the  core-halo  simulation  beam.  The  simulation  reproduces  the 
n,  rms 
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radial  oscillation  frequency  and  the  damping  rate  observed  in  the 
experiments.  Total  emittance  growth  in  the  simulation  was  29%. 

Experiments  were  also  done  using  the  three-element  match  of  Ref.  21. 
Experimental  major  and  minor  radii,  again  defined  in  terms  of  the  time- 
averaged  current-density  contour  at  10%  of  the  peak  current  density,  are 
plotted  versus  transport  distance  in  Fig.  23.  In  this  case,  the  mismatch 
oscillations  are  greatly  reduced. 

Simulation  results  using  the  core-halo  distribution  as  in  Fig.  22  are 
shown  in  Fig.  24.  Here  we  have  used  a  single-element  match  with  peak 
quadrupole  gradient  b^  =  75  G/cm  located  5  cm  from  the  stellarator 
entrance.  As  before,  the  axial  lens  profile  is  that  shown  in  Fig.  15.  We 
find  smaller  radial  oscillations,  as  in  the  experiment,  and  emittance 
growth  of  only  11%. 


24 


VI.  CONCLUSIONS 


Simulation  studies  of  recent  SLIA  experiments  and  of  the  SLIA  POCE 
have  been  performed.  Simulation  results  shoved  that  as  the  intense  beam  is 
transported  through  each  bend,  the  beam  experiences  self  electro-  and 
magnetostatic  forces.  To  account  for  these  forces,  the  design  value  of  the 
vertical  bending  field  had  to  be  increased  by  as  much  as  25%.  Optimum 
vertical  field  values,  which  were  in  reasonable  agreement  with  the 
theoretical  analysis  of  Ref.  6,  were  determined  for  each  bend. 
Furthermore,  we  found  that  this  correction  to  the  vertical  field  did  not 
adversely  effect  the  chromatic  properties  of  the  bend. 

Simulations  of  transport  into  the  stellarator  fields  shoved  that 
radial  mismatch  oscillations  are  likely  to  be  excited  both  at  the 
stellarator  entrance  and  at  the  exit.  Simulation  results  were  quite 
sensitive  to  the  spatial  distribution  of  the  beam  particles.  Ve  found 
that,  if  the  radial  profile  of  the  beam  is  rounded,  the  mismatch 
oscillations  damp  and  that  this  damping  is  accompanied  by  emittance  growth. 
These  phenomena  did  not  occur  when  the  beam  was  initialized  with  a 
Kapchinskij-Vladimirskij  distribution,  however.  Ve  believe  that  the 
observed  damping  of  radial  mismatch  oscillations  is  due  to  phase  mixing 
caused  by  nonlinear  space-charge  forces  associated  with  a  non-uniform  beam 
profile. [14-16]  By  contrast,  the  K-V  distribution  has  a  uniform  spatial 
profile. 

We  have  presented  a  novel  method  of  matching  the  beam  onto  and  off  of 
the  stellarator  fields.  Here  we  use  a  single  thick  quadrupole  lens  on  each 
end  of  the  stellarator  to  transform  the  circular  cross-section  of  the  beam 
in  the  solenoid  field  into  an  elliptical  shape  that  closely  approximates 
the  equilibrium  configuration  in  the  stellarator  fields.  The  lens  strength 
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and  orientation  are  adjusted  to  obtain  the  correct  eccentricity  and  phase 
at  the  stellarator  entrance.  An  identical  lens,  with  only  its  orientation 
adjusted,  is  used  at  the  exit  to  bring  the  beam  back  to  a  circular  cross- 
section.  Simulations  show  that  matching  can  be  accomplished  in  this  way 
and  that  the  degree  of  emittance  growth  can  be  greatly  reduced  by  the  use 
of  such  single-element  matching  lenses. 

We  also  performed  end-to-end  simulations  of  the  SLIA  POCE,  where  the 
beam  was  taken  through  all  three  bends.  These  showed  that  matching 
elements  can  reduce  total  emittance  growth  from  52 X  to  10%.  While  the 
level  of  emittance  growth  found  in  the  end-to-end  simulations  was  modest 
even  when  no  matching  elements  were  used,  such  elements  will  certainly  be 
necessary  in  a  50-100  MeV  SLIA  design,  where  many  more  bends  would  be 
required. 

We  have  also  explored  other  approaches  to  the  SLIA  matching  and 
transport  problems,  which  are  not  reported  here.  For  instance,  discrete, 
alternating-gradient  quadrupole  lenses  were  considered  as  an  alternative  to 
stellarator  windings.  The  motivation  for  this  study  was  the  possibility  of 
improved  stability  with  regard  to  the  electromagnetic  instability  of  Refs. 
30-32.  It  was  found,  however,  that  alternating-gradient  transport  is 
comparably  unstable. [33,34]  It  was  further  found  that  the  instability 
could  be  avoided  in  the  stellarator  configuration  by  carefully  choosing 
system  parameters. [31,32,35]  We  also  considered  a  design  in  which 
stellarator  windings  were  used  on  the  straight  sections  as  well  as  on  the 
bends.  This  design  does  not  avoid  the  complication  of  changing  the 
parameters  of  the  stellarator  and  solenoid  fields  between  each  bend,  in 
order  to  keep  the  bends  achromatic  as  the  beam  energy  increases.  We  found 
that  the  stellarator-stellarator  transitions  (i.e.,  between  the  parameters 
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of  bends  one  and  two)  were  no  less  problematical  than  the  solenoid- 
stellarator  transitions  of  the  POCE.  This  difficulty  persisted  even  for 
designs  in  which  each  bend  had  the  same  stellarator  pitch  length. 

Finally,  simulations  were  performed  to  model  and  understand  the  recent 
SLIA  transport  experiments. [8-10]  A  core-halo  distribution,  similar  to 
that  of  the  850  keV,  200  A  experimental  beam,  was  used  in  simulations  of 
transport  through  the  stellarator  fields,  with  and  without  matching  lenses. 
These  simulations  reproduced  each  of  the  key  features  of  the  experiments, 
including  the  motion  of  the  beam  cross-section,  which  followed  the  pitch  of 
the  stellarator  windings  in  all  cases,  the  amplitude,  frequency  and  damping 
rate  of  the  radial  mismatch  oscillations,  and  the  reduction  in  the 
amplitude  of  these  oscillations  through  the  use  of  matching  elements  at  the 
stellarator  entrance  and  exit. 

In  the  850  keV  simulations,  the  sensitivity  of  the  results  to  the 
spatial  distribution  of  the  beam  was  of  key  importance  in  our  efforts  to 
reproduce  the  experimental  results.  Specifically,  the  initial  measurements 
of  the  spatial  profile  of  the  beam  showed  only  the  central  core. 
Simulation  results  with  the  same  profile  (Fig.  19)  showed  a  25X  disparity 
in  the  frequency  of  the  mismatch  oscillations  between  simulation  and 
experiment  (see  Fig.  18).  By  adding  a  halo,  we  found  that  this  disparity 
could  be  eliminated.  The  addition  of  a  halo  was  motivated  by  the  fact  that 
such  core-halo  distributions  are  associated  with  the  damping  of  radial 
mismatch  oscillations[29]  as  might  be  present  in  the  injector. [10,201  The 
existence  of  the  halo  was  later  confirmed  experimentally. [ 10] 
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APPENDIX:  FIRST  ORDER  FIELD  CORRECTIONS  IN  A  CURVED  GEOMETRY 

In  the  ELBA[11]  and  SST[ 13 ]  particle  simulation  codes,  the 
electromagnetic  fields  are  computed  in  terms  of  the  scalar  and  vector 
potentials  using  cylindrical  coordinates.  In  all  cases,  the  beam  is 
assumed  to  be  propagating  within  a  conducting  pipe  of  radius  rv-  Following 
the  approach  originated  by  Lee[36],  the  equations  are  solved  using  the 
gauge 


V*L  *  0  »  (A1) 

where  A^  is  the  component  of  the  vector  potential  perpendicular  to  the 
direction  of  propagation.  This  gauge  works  well  for  problems  in  which  the 
beam  is  nearly  axisymmetric  and  dominates  the  field  production.  In  the 
ELBA  code,  a  change  of  variables  is  made  from  (r,0,z,t)  to  (r,9, C*ct- 
z,T«t).  The  (C, t)  variables  are  numerically  advantageous  in  that  a)  a 
relativistic  beam  is  nearly  stationary  in  these  coordinates  and  b) 
information  flow  is  always  towards  larger  values  of  C-  The  SST  code  makes 
the  additional  simplifying  assumption  that  9/9C  «  0.  This  approximation  is 
useful  when  self-forces  arising  from  C-derivatives  of  the  beam  density  and 
current  can  be  neglected. 

In  terms  of  the  potentials,  A,,  *  and  a  ■  A  -  Maxwell's  equations 
in  the  gauge  of  Eq.  (Al)  can  be  written  as [11, 37] 

3E 

vj*  ♦  -4np  ,  (A2) 

3E 

V^a  +  -  -4n(Jz/c  -  p)  (A3) 
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and 


4n  T  n  (3a  1  3^| 

-~il  -  'ite'  c  F?J  ’ 


where 


E 

z 


1  3* 
c  3t  * 


(A4) 


(A5) 


For  propagation  around  a  bend,  ve  neglect  corrections  to  A  and  A.  and 

V  v 

retain  terms  in  the  equations  for  a  and  *  to  first  order  in  r  /R,  where  R 

w 

is  the  radius  of  curvature  of  the  bend.  The  corrections  are  treated  as 
explicit  terms  in  the  solution  algorithms.  Maxwell's  equations  for  the 
potentials  in  the  bend  are  derived  by  transforming  to  a  local  orthogonal 
coordinate  system  (r,0,s)  which  is  centered  on  the  axis  of  the  beam  pipe 
with  s  measuring  arc  length  along  the  bend.  The  (r,0,s)  coordinates  can  be 
written  in  terms  of  global  Cartesian  coordinates  (x,y,z)  with  the  origin  at 
the  center  of  the  bend: 


x  «  (R  +  r  cos0)cos(s/R)  , 


(A6a) 


y  -  r  sin0 


(A6b) 


and 


z  *  (R  +  r  cos0)sin(s/R)  ,  (A6c) 

for  r  <  r  . 

-  v 

Maxwell's  equations  for  a  and  *  are  written  in  the  (r,0,s)  coordinate 
system  and  expanded  to  first  order  in  r/R.  We  then  make  a  change  of 

variables  from  (r,0,s,t)  to  (r , 0, C*ct-s , ) .  The  result  is 
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where 


+  jz  “  -  Hcose  If  "  H +  c cose  b&z 


a 2  a2  n  *K\ 

c  cose  FEK  +  r  COS0  Jjf  +  COs9(f  FE  +  ) 
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ii 

3C 


(A9) 


and  Jg  is  the  current  density  in  the  s  direction.  Because  we  neglect 

corrections  to  A  and  Aa,  the  gauge  condition,  Eq.  (Al),  and  Eq.  (A4) 
r  u 

remain  unchanged.  The  remaining  electric  and  magnetic  field  components  are 
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Table  I. 


Parameters  for  the  SLIA  POCE:  bend  design  energy  E,  axial  field 
Bg,  quadrupole  gradient  b^  and  quadrupole  pitch  length  X^. 


Bend 

E  (MeV) 

Bs  (G) 

bQ  (G/cm) 

Xp  (cm) 

1 

4.0 

5492. 

251.6 

68  33 

2 

5.5 

5335. 

207.0 

83.78 

3 

7.0 

5483. 

182.3 

125.67 

Table  II.  Numerically  determined  optimum  vertical  field  values  for  each 
POCE  bend.  Values  are  in  units  of  Gauss.  For  comparison,  the 
ideal  values  (no  hoop  forces)  and  the  predicted  values  from  Eq. 
(2)  are  also  given. 


Bend  By 


Bv  [ideal]  Bv  [4-1]  By  [4-2] 


1  233.0  186.9 

2  296.0  249.5 

3  354.5  312.5 


236.2  248.6 
^98.6  311.0 
.'-0.1  372.5 


Table  III.  Location  for  each  bend,  each  ideal  acceleration  and  each 
solenoid  transition  for  the  end-to-end  POCE  simulations. 

Lengths  are  in  centimeters.  Note  that  the  acceleration  after 
the  third  bend  is  not  included  and  that  the  POCE  design  features 
longer  straight  sections  than  those  used  here. 

Bend  Bend  location  Acceleration  Solenoid  transition 


1 

2 

3 


50  -  301.3  400  -  450 
550  -  801.3  900  -  950 
1050  -  1301.3 


422.5  -  427.5 

922.5  -  927.5 


37 


xeai  o 


0.25  : - " - 

0.20 : 

- 

&  °-15{  n 

*  ) 

s  0.10L  □ 

[ 

I 

0.05  ^ 

h 

L 

L 

0  . _ 1. _ . _ _ 

-10  -5 


Fig.  7. 


Beam  centroid  displacement 
error  in  the  axial  magneti 


mismatch  AE  *  -102. 


0  5  10 

ABS  (%) 

at  the  bend  exit  Sr  versus  percentage 

field  AB  for  runs  with  beam  energy 
s 


44 


cm 


0.10  h 


h 

i 

L 


f- 


0.05  r 


Abq  (%) 


Fig.  8.  Beam  centroid  displacement  at  the  bend  exit  &r  versus  percentage 
error  in  the  stellarator  gradient  dbq  for  runs  with  beam  energy 
mismatch  AE  -  -10%. 
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Major  and  minor  radii,  r+  (solid)  and  r  (dashed),  of  the  beam 
cross-section  versus  transport  distance  for  beams  with  various 
profiles:  (a)  K-V,  (b)  parabolic,  (c)  Gaussian  and  (d)  Bennett. 
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Fig.  10. 
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Emictance  t  _  versus  transport  distance  for  each  of  the  four 
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cases  of  Fig.  9:  K-V  (solid  line),  parabolic  (dashed),  Gaussian 

(dotted)  and  Bennett  (long  dashes). 
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Fig.  11.  Envelope  code  results  showing  major  and  minor  beam  radii  versus 
transport  distance  for  several  quadrupole  lens  lengths  L^:  a)  LQ 
-  1  cm,  b)  Lq  .  2.5  cm,  c)  LQ  «  5  cm  and  d)  LQ  *  10  cm.  Here, 
LQbQ  »  2000  G  in  all  cases,  where  bQ  is  the  peak  quadrupole 
gradient  of  the  lens.  Beam  parameters  correspond  to  those  of  the 
first  P0CE  bend. 
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Fig.  12.  Envelope  code  results  shoving  major  and  minor  beam  radii  versus 
transport  distance  for  a  quadrupole  lens  with  length  LQ  *  5  cm 
and  quadrupole  gradient  bg  «  400  G/cm  as  in  Fig  11  (c).  Beam 

parameters  correspond  to  those  of  the  third  POCE  bend. 
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Fig.  13.  Simulation  result  shoving  major  and  minor  beam  radii,  r+  and  r  , 
as  a  function  of  transport  distance.  Here,  single-element 
matching  has  been  used. 
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Fig.  14.  Simulation  result  showing  normalized  RMS  emittance  as  a  function 
of  transport  distance  for  the  single-element  matching  case  as 
shown  in  Fig.  13  (solid  line)  and  for  the  case  with  no  matching 
as  shown  in  Fig.  9(b)  (dashed  line). 
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Fig.  17.  End-to-end  POCE  simulation  results  showing  (a)  major  and  minor 
radii  and  (b)  emittance  versus  transport  distance.  Single¬ 
element  matching  lenses  are  used  at  the  entrance  and  exit  of  each 
bend.  The  overall  emittance  growth  is  10%. 
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Fig.  20.  Experimental  data  from  Tiefenback  et  al.[10]  showing  charge 
density  versus  transverse  coordinate  x. 
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Fig.  22.  Simulation  result  shoving  (a)  major  and  minor  beam  radii  and  (b) 
normalized  RMS  emittance  as  a  function  of  transport  distance  for 
a  beam  with  a  core-halo  distribution. 

59 


n ,  r  i  n  s 


0.5  ~ 


0  L_ _ : . — _ * _ . . - _ . _ 

0  100  200  300  400  500 

z  (cm) 

Fig.  24.  Simulation  result  shoving  (a)  major  and  minor  beam  radii  and  (b) 
normalized  RMS  emittance  as  a  function  of  transport  distance.  In 
this  case,  a  single-element  matching  lens  was  used. 
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